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We numerically construct asymptotically anti-de Sitter (AdS) black holes in four dimensions that 
contain only a single Killing vector field. These solutions, which we coin black resonators, link 
the superradiant instability of Kerr-AdS to the nonlinear weakly turbulent instability of AdS by 
connecting the onset of the superradiance instability to smooth, horizonless geometries called geons. 
Furthermore, they demonstrate non-uniqueness of Kerr-AdS by sharing asymptotic charges. Where 
black resonators coexist with Kerr-AdS, we find that the black resonators have higher entropy. 
Nevertheless, we show that black resonators are unstable and comment on the implications for the 
endpoint of the superradiant instability. 


Introduction. As the simplest of gravitating objects, 
black holes (BHs) play a fundamental role in our under¬ 
standing of general relativity. Indeed, four-dimensional, 
asymptotically flat BHs are stable and uniquely speci¬ 
fied by their asymptotic charges However, there are 
circumstances where stability and uniqueness can be vio¬ 
lated, such as those in higher dimensions HE]. We will 
argue that this can also be accomplished in four dimen¬ 
sions with asymptotically anti-de Sitter (AdS) BHs. 

Unlike Minkowski or de Sitter space, AdS contains a 
timelike boundary at conformal infinity where reflect¬ 
ing (energy and angular momentum conserving) bound¬ 
ary conditions are typically imposed to render the ini¬ 
tial value problem well posed [12]. The presence of this 
boundary has drastic consequences for the stability of so¬ 
lutions in AdS. For example, rotating BHs contain an er- 
goregion from which energy can be extracted by the Pen¬ 
rose process [IJj . For waves, this phenomenon is called 
superradiance [IdUIH] . In AdS, these waves return after 
scattering from the boundary and extract more energy. 
The process continues until the waves contain enough en¬ 
ergy to backreact on the geometry, causing the so-called 
superradiant instability [nHisi- 

The reflecting boundary also has implications for the 
stability of AdS itself. A nonlinear instability may occur 
if an excitation with arbitrarily small, but finite energy 
around AdS continues to reflect off the boundary and 
eventually forms a BH. There is numerical evidence in 
support of this instability with a spherically symmetric 
scalar field . There is additionally a proposed per¬ 

turbative explanation for this instability |20] which ap¬ 
plies to pure gravity and beyond spherical symmetry |28j . 
At linear order in perturbation theory, AdS contains an 
infinite tower of evenly-spaced normal modes. At higher 
orders, resonances between modes cause higher modes to 
be excited that grow linearly in time. In the generic case, 
this leads to a breakdown of perturbation theory, and is 
interpreted as the beginnings of a nonlinear instability. 
This instability is called weakly turbulent due to this en¬ 


ergy shift from longer to shorter length scales. 

Though there is a breakdown of perturbation theory 
for generic initial data, perturbation theory survives to 
arbitrarily high orders when only a single mode is excited. 
This leads to a family of horizonless time-periodic solu¬ 
tions called oscillons {boson stars) for a real (complex) 
scalar field [201 ED and geons for pure gravity [281124j . 
These geons can be thought of as nonlinear normal modes 
of AdS and are solutions that contain only a single Killing 
field. Any gravitational radiation emitted by the geon is 
balanced by absorption of waves reflected from the AdS 
boundary. Since perturbation theory breaks down for 
two-mode initial data, geons can be thought of as a basis 
for the nonlinear instability [2S1 [25] . 

We will construct a new family of BHs that joins the 
onset of the superradiant instability of Kerr-AdS [mui] 
[40] to the geons. Since these BHs are time-periodic and 
single out a particular frequency, we call them black res¬ 
onators. One limit of black resonators corresponds to the 
onset of the superradiant instability. Black resonators 
are thus the BHs predicted by [2S] and alluded to in 
[IIl|23l|26]. The opposite, zero-size limit of black res¬ 
onators corresponds to the geons where this frequency is 
given by a nonlinear normal mode of AdS [22] [25] . Like 
the geons, black resonators are time-periodic and have 
only one Killing field. The Killing field is also the hori¬ 
zon generator, so black resonators have bifurcate Killing 
horizons and globally well-defined horizon temperatures. 

Numerical Approach. We will search for solutions 
to the Einstein equation with a negative cosmological 
constant A = —Zjl?. That is, we will solve 
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where L is the AdS length scale, gab the metric and Rab 
the Ricci tensor. We wish to find asymptotically global 
AdS solutions with a boundary metric conformal to the 
Einstein static universe 

ds^dy = —dt^ -I- d9^ sin^ 9dip^ . (2) 
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We use the DeTurck method m which proceeds by 
writing down any reference metric g of our choice that 
shares the same causal structure as the solution g we wish 
to find. Then, rather than solve the Einstein equation 
0 . we instead solve the Einstein-DeTurck equation: 

Rab + ^ ~ ^ (aCh) = 0 i (3) 

with where r(g) is the Christof- 

fel symbol associated with a metric g. Since we are solv¬ 
ing different equations, we verify a posteriori that our 
solutions to ([^ satisfies = 0 to machine precision, and 
hence is also a solution to Q. In some cases, it is pos¬ 
sible to prove that all solutions of ^ must have = 0 
[25], but we do not have such a proof for our case of in¬ 
terest. However, we have verified that ([^ yields elliptic 
partial differential equations for which local uniqneness 
theorems exist |2^. Thus, solutions with ^ 7 ^ 0 cannot 


have ^ arbitrarily small for all ranges of parameters, and 
are therefore distinguishable from solutions with ^ = 0 . 

Following the DeTurck method, we need an ansatz and 
a reference metric containing a single Killing field, 

K = dt + , (4) 

where fin will be the horizon angular velocity. KT is a 
Killing field, but dt and are not. Thus, the solution 
will neither be time independent nor axisymmetric, but 
is instead time-periodic. We now perform the following 
change of variables: dr = dt and dcj) = dip + ilndt. In 
these coordinates, K = dr and the boundary metric is: 

ds^dy = —dr^ -I- dO^ + sin^ 9{d(j) — ^l^dr)^ . (5) 

For a general ansatz containing a Killing horizon gen¬ 
erated by K, we choose: 


ds^ = 




(l-y2)2 


- y^qiMy) (dr -h y qedy) 
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(dr -f yxV'2 - grdy + y^x^/2 - x^ ggdr)' 


rVS — x'^qgdx 


+ (1 - xMy+yi ( M - y^lsdr -h ^ + y qiody 


( 6 ) 


with A(y) = (I - vM + y+(3 - -f y'^) and qi = 
qi{x,y,(j)). As a reference metric, we take 0 with 
91 = 92 = 93 = 94 = 1 , 95 = and 9 * = 0 , for i > 6 . 

If D// = 0, our reference metric is the Schwarzschild- 
AdS BH with entropy S = ny^L"^ and temperature 


l + 3y^ 

47ry+L 


(7) 


To recover more familiar coordinates, redefine rjL — 
y+/(l — y^) and cos0 = X'<j2 — r^. We chose this ref¬ 
erence metric rather than Kerr-AdS for tidiness. 

Since black resonators branch from the onset of the 
superradiant instability of Kerr-AdS nms], let us now 
describe this instability in more detail. Perturbations of 
Kerr-AdS are labeled by a type (scalar or vector) and 
polar and azimutal wavenumbers ^ and m, respectively. 
While all Kerr-AdS BHs with VIhL > 1 are superradi¬ 
ant unstable, different perturbative modes have onsets in 
difference places in parameter space. At these onsets are 
zero modes which mark the appearance of a new family 
of solutions which will be the black resonators. 

For definitiveness and simplicity, we will focus on the 
scalar mode £ = m = 2, which has an extra discrete 
symmetry under a; —> —x. The Kerr-AdS parameters 
corresponding to the onset of the superradiant instability 


of this mode were obtained in [nidg] m). 

Let us now discuss boundary conditions. At y = 0, we 
demand a regular bifurcate Killing horizon generated by 
dr with temperature 0 . This amounts to 9 i(x, 0,(/)) = 
92 ( 2 :, 0, 4>), with the remaining functions having Neumann 
conditions. The boundary (y = 1) must approach that of 
the reference metric, which is a Dirichlet condition. The 
discrete symmetry x —>■ —x requires that all of the qt’s 
have Neumann boundary conditions at a: = 0. At a; = 1, 
regularity requires qT(\,y,(j)) = q^(\,y,(p) = 0 and Neu¬ 
mann conditions for the remaining functions. Finally, 
m = 2 requires that (j) be periodic in 0 G ( 0 , tt]. 

With equations and boundary conditions, we can 
now solve the system of partial differential equations 
by numerical methods. We use a standard Newton- 
Raphson algorithm and discretise the Einstein-DeTurck 
equations using pseudospectral collocation (Chebyshev- 
Gauss-Lobatto nodes along the x and y directions, and 
Fourier nodes along the (j) direction). The resulting alge¬ 
braic linear systems are solved by LU decomposition. 

Our reference metric is not Kerr-AdS, so we must con¬ 
struct Kerr-AdS numerically. We fix y+, assume inde¬ 
pendence in (j), and slowly increase VIh from zero until we 
are at the onset of the superradiant instability. Note that 
this is consistent with our boundary conditions above. 
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FIG. 1: E vs. J phase diagram of AdS black holes with Kerr- 
AdS BHs (all regions above thick solid black line), geons in 
|24 | (bottom yellow dashed line), and black resonators (data 
points). Also plotted are Kerr-AdS BHs with Hr/L < 1 
(dashed purple line and purple region above it), extremal 
Kerr-AdS BHs (thick solid black line), and the onset of scalar 
m = i = 2 mode in |19| (thin blue line). 


From here, we have tried to obtain black resonators 
by varying Uh and using a perturbed Kerr-AdS as 
a seed. This proved unsuccessful since Kerr-AdS is too 
strong an attractor for Newton-Raphson. Instead, we 
used y+ and a wiggliness parameter e defined as 

pTT 

£= < 79 ( 1 , 0 , (/>) sin(m (/))d(/). (8) 

Jo 

We add the defining equation Q to our system of equa¬ 
tions and fin as an extra unknown. Any solution with 
e ^ 0 has (j) dependence and hence cannot be Kerr-AdS. 

Results. Our main result is shown in Fig.where the 
energy [E) versus angular momentum (J) phase diagram 
of rotating BHs in AdS is presented in units of the AdS 
radius L. Kerr-AdS BHs lie above the thick solid black 
line, which refers to extremal (T = 0) Kerr-AdS BHs. 
Kerr-AdS BHs with < I are likely linearly stable 
m ISO] and lie in upper-left region above the dashed 
purple line (the purple line itself refers to flnL = 1). 
The onset of the superradiant instability for the scalar 
i = m = 2 modes [19] lie on the thin blue line. The 
geons from [21] lie on the bottom dashed line. The data 
points refer to the black resonators we have constructed, 
with parameters indicated by the plot markers. 

From Fig.j^ we see that our solutions at small e recover 
the onset curve of m- This is a strong consistency check 
since the onset curve was generated using the Teukolsky 
equation in AdS, which only indirectly uses the metric 
[18| . All of our black resonators also have flnL > 1, 
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FIG. 2: Gomparison of numerical e = 0.1 black resonators 
(blue disks) and the pertnrbative result Q (red diamonds). 


which indicates that no Killing vector can be found near 
the conformal boundary that is everywhere timelike. We 
also see that some black resonators extend below the ex¬ 
tremal limit of Kerr-AdS BHs and are therefore the only 
known regular solutions with these asymptotic charges. 

Though we focused on the scalar £ = m = 2 mode, a 
perturbative expansion can predict the entropy of small 
E and J black resonators for more general modes m- 
For scalar modes with £ = m, the entropy is 


S = AttE^ 




—r 

eI, 


(9) 


For small black resonators, Fig.j^shows good agreement 
between our numerical results and the perturbative pre¬ 
diction of ([^ with m = 2. Since ([^ assumes that zero- 
size black resonators merge with geons, this agreement 
can be taken as evidence that black resonators connect to 
geons. Fig. shows black resonators approaching geons, 
but only down to size y+ ~ 0.07. 

We note that if black resonators with arbitrary m are 
connected to geons, the minimum E for fixed J that black 
resonators can have occurs when m —>■ -boo for an arbi¬ 
trarily small black resonator {i.e. a geon with m = -boo). 
This configuration saturates the BPS bound E = JfL. 

Another quantity of interest is the entropy of black res¬ 
onators compared with that of Kerr-AdS. At the same 
asymptotic charges E and J, we find that black res¬ 
onators have higher entropy than Kerr-AdS ( see for in¬ 
stance Fig. [^. These solutions merge at the onset of 
superradiance through a second order phase transition. 

We show the energy density on the 5^ of the boundary 
metric (§ in Fig. 1^ [31]. This figure represents an instant 
in time and should be imagined as rotating with angular 
velocity Hh (*-e. in a time-periodic way). 

Outlook. We have constructed new BHs in AdS with 
a single Killing field which we call black resonators [42j . 
Black resonators are not ruled out by Hawking’s rigidity 
theorem [5M55] since the single Killing held also gener¬ 
ates the horizon. These BHs branch from the onset of the 
superradiant instability in Kerr-AdS, and extend in their 
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FIG. 3: Entropy of rotating BHs in AdS versus their energy 
with Kerr-AdS BHs (solid blue line) and black resonators (red 
triangles). The black resonators all have y+ — 0.16. 
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FIG. 4: Energy density in units of L at a moment in time of 
a black resonator with y+ = 0.16 and e = 0.1. 




zero-size limit to smooth horizonless geometries called 
geons [23l |24j . The existence of black resonators proves 
that Kerr-AdS is non-unique, even in four dimensions. 
We focused on the scalar progenitor mode m = £ = 2, 
but expect similar behaviour for other m and i. This 
would mean a countably infinite violation of uniqueness 
for rotating BHs in AdS with flnL > 1 and E > J/L. 

In retrospect, new BHs could have been anticipated 
from the AdS/CFT correspondence. Since CFTs are ex¬ 
pected to saturate the bound E > J/L, but Kerr-AdS 
BHs do not, another BH might fill the gap. Though, this 
argument does not suggest that these BHs have a single 
Killing field or are connected to the superradiant onset. 

The precise boundary CFT interpretation of these in¬ 
stabilities and the black resonators remain mysterious. 
We note that superradiance is not particular to four di¬ 
mensions, and occurs also in AdS^. Furthermore, includ¬ 
ing the full AdS^ x bulk geometry, so that the bound¬ 
ary field theory is specifically N = A super Yang-Mills, 
does not cure this instability. 

Though these black resonators have more entropy than 
Kerr-AdS, we argue that they are unstable, so they can¬ 


not be the endpoint of the superradiant instability. While 
black resonators with progenitor modes £ = m = 2 should 
be stable to perturbations with m = £ = 2, they are un¬ 
stable to higher m modes. The reason for this is that 
small black resonators are well-approximated (as con¬ 
firmed in ([^ and Fi g. [2 ) by a small Kerr-AdS BH at 
the centre of a geon mm, and small Kerr-AdS BHs 
are still unstable to higher m modes. More precisely, the 
results of |36j mathematically prove that our solutions are 
unstable, since no Killing vector held that is everywhere 
timelike can be found at the conformal boundary. 

Our results support the conjecture of m that there is 
no stationary endpoint to the superradiant instability in 
AdS. Instead, modes with increasing m continue to be ex¬ 
cited and develop. It may be possible for additional small 
BHs to form as energy is deposited into higher m modes. 
While such conhgurations should exist, they are them¬ 
selves superradiantly unstable [38l [39]. Though classical 
evolution may continue indehnitely, eventually the in¬ 
creasingly high m modes will reach sub-Planckian length 
scales. This may be viewed as a violation of the spirit of 
cosmic censorship in that initial data well-described clas¬ 
sically leads to a situation requiring quantum mechanics. 

Acknowledgments 

It is a pleasure to thank Dave Berenstein and Don 
Marolf for discussions. The authors thankfully acknowl¬ 
edge the computer resources, technical expertise, and 
assistance provided by CENTRA/IST. Some compu¬ 
tations were performed at the cluster ‘Baltasar-Sete- 
Sois’, supported by the DyBHo-256667 ERG Starting 
Grant; and on the GOSMOS Shared Memory system at 
DAMTP, University of Gambridge operated on behalf of 
the STFG DiRAG HPG Facility and funded by BIS Na¬ 
tional E-infrastructure capital grant ST/J005673/I and 
STFG grants ST/H008586/1, ST/K00333X/1. O.J.G.D. 
is supported by the STFG Ernest Rutherford grants 
ST/K005391/I and ST/M004147/1. This research re¬ 
ceived funding from the European Research Gouncil un¬ 
der the European Gommunity’s 7th Framework Pro¬ 
gramme (FP7/2007-2013)/ERG grant agreement no. 
[247252]. B.W. is supported by European Research 
Gouncil grant no. ERG-2011-StG 279363-HiDGR. 


* Electronic address: 0.J.Gampos-Dias@soton.ac.uk 
^ Electronic address; J.E.Santos@damtp.cam.ac.uk 

* Electronic address; B.Way@damtp.cam.ac.uk 

[1] D. Robinson, Four decades of black holes uniqueness 
theorems, in The Kerr Spacetime: Rotating Black Holes 
in General Relativity (D. L. Wiltshire, M. Visser, S. M. 
Scott, ed.). Gambridge University Press, 2009. 

[2] R. Gregory and R. Laflamme, Black strings and 
p-branes are unstable, Phys.Rev.Lett. 70 (1993) 
2837-2840, hep-th/9301052 . 










5 


[3] L. Lehner and F. Pretorius, Black Strings, Low 
Viscosity Fluids, and Violation of Cosmic Censorship, 
Phys.Rev.Lett. 105 (2010) 101102, arXiv: 1006.5960 . 

[4] R. Emparan and H. S. Reall, A Rotating black ring 
solution in five-dimensions, Phys.Rev.Lett. 88 (2002) 
101101, hep-th/0110260 . 

[5] O. J. Dias, P. Figueras, R. Monteiro, J. E. Santos, and 
R. Emparan, Instability and new phases of 
higher-dimensional rotating black holes, Phys.Rev. D80 
(2009) 111701, arXiv:0907.2248 . 

[6] M. Shibata and H. Yoshino, Nonaxisymmetric instability 
of rapidly rotating black hole in five dimensions, 

Phys.Rev. D81 (2010) 021501, arXiv:0912.3606. 

[7] M. Shibata and H. Yoshino, Bar-mode instability of 
rapidly spinning black hole in higher dimensions: 
Numerical simulation in general relativity, Phys. Rev. 
D81 (2010) 104035, arXiv: 1004.4970 . 

[8] O. J. Dias, P. Figueras, R. Monteiro, and J. E. Santos, 
Ultraspinning instability of rotating black holes, 

Phys.Rev. D82 (2010) 104025, arXiv: 1006.1904 . 

[9] O. J. Dias, P. Figueras, R. Monteiro, and J. E. Santos, 
Ultraspinning instability of anti-de Sitter black holes, 
JHEP 1012 (2010) 067, arXiv: 1011.0996 . 

[10] O. J. Dias, J. E. Santos, and B. Way, Rings, Ripples, 
and Rotation: Connecting Black Holes to Black Rings, 
JHEP 1407 (2014) 045, arXiv: 1402.6345 . 

[11] J. E. Santos and B. Way, The Black Ring is Unstable, 
arXiv:1503.0072 

[12] H. Friedrich, Einstein equations and conformal structure 
- Existence of anti de Sitter type space-times, 
J.Geom.Phys. 17 (1995) 125-184. 

[13] R. Penrose and R. Floyd, Extraction of rotational 
energy from a black hole. Nature 229 (1971) 177-179. 

[14] Y. B. Zeldovich, Generation of Waves by a Rotating 
Body, JETP Lett. 14 (1971) 180. 

[15] A. A. Starobinsky and S. M. Churilov, Amplification of 
electromagnetic and gravitational waves scattered by a 
rotating black hole, Sov. Phys. JETP 38 (1973) 1. 

[16] S. Teukolsky and W. Press, Perturbations of a rotating 
black hole. Ill - Interaction of the hole with 
gravitational and electromagnetic radiation, 

Astrophys.J. 193 (1974) 443-461. 

[17] V. Cardoso and O. J. Dias, Small Kerr-anti-de Sitter 
black holes are unstable, Phys.Rev. D70 (2004) 084011, 
hep-th/0405006 . 

[18] Dias, Oscar J.C. and Santos, Jorge E., Boundary 
Conditions for Kerr-AdS Perturbations, JHEP 1310 
(2013) 156, arXiv: 1302.1580. 

[19] V. Cardoso, O. J. Dias, G. S. Hartnett, L. Lehner, and 
J. E. Santos, Holographic thermalization, quasinormal 
modes and superradiance in Kerr-AdS, JHEP 1404 
(2014) 183, arXiv:1312.5323 . 

[20] P. Bizon and A. Rostworowski, On weakly turbulent 
instability of anti-de Sitter space, Phys.Rev.Lett. 107 
(2011) 031102, arXiv:1104.3702 . 

[21] A. Buchel, L. Lehner, and S. L. Liebling, Scalar 
Collapse in AdS, Phys.Rev. D86 (2012) 123011, 
arXiv:1210.0890 . 

[22] H. Okawa, J. C. Lopes, and V. Cardoso, Collapse of 
massive fields in anti-de Sitter spacetime, 

arXiv:1504.0520 

[23] O. J. Dias, G. T. Horowitz, and J. E. Santos, 
Gravitational Turbulent Instability of Anti-de Sitter 
Space, Class. Quant. Grav. 29 (2012) 194002, 


arXiv: 1109.1825. 

[24] G. T. Horowitz and J. E. Santos, Geons and the 
Instability of Anti-de Sitter Spacetime, 

arXiv:1408.5906 

[25] H. S. Reall, Higher dimensional black holes and 
supersymmetry, Phys.Rev. D68 (2003) 024024, 
hep-th/0211290. 

[26] H. K. Kunduri, J. Lucietti, and H. S. Reall, 
Gravitational perturbations of higher dimensional 
rotating black holes: Tensor perturbations, Phys. Rev. 
D74 (2006) 084021, hep-th/0606076 . 

[27] M. Headrick, S. Kitchen, and T. Wiseman, A New 
approach to static numerical relativity, and its 
application to Kaluza-Klein black holes. 

Class. Quant. Grav. 27 (2010) 035002, 

[arXiv:0905.1822. 

[28] P. Figueras, J. Lucietti, and T. Wiseman, Ricci solitons, 
Ricci flow, and strongly coupled GET in the 
Schwarzschild Unruh or Boulware vacua. 

Class. Quant. Grav. 28 (2011) 215018, 
arXiv:1104.4489 . 

[29] A. Adam, S. Kitchen, and T. Wiseman, A numerical 
approach to finding general stationary vacuum black 
holes. Class. Quant. Grav. 29 (2012) 165002, 
arXiv:1105.6347 . 

[30] S. Hawking and H. Reall, Charged and rotating AdS 
black holes and their CRT duals, Phys.Rev. D61 (2000) 
024014, hep-th/9908109 . 

[31] S. de Haro, S. N. Solodukhin, and K. Skenderis, 
Holographic reconstruction of space-time and 
renormalization in the AdS / CFT correspondence, 
Commun.Math.Phys. 217 (2001) 595-622, 
hep-th/0002230. 

[32] S. Hawking, Black holes in general relativity, 
Commun.Math.Phys. 25 (1972) 152-166. 

[33] S. Hollands, A. Ishibashi, and R. M. Wald, A Higher 
dimensional stationary rotating black hole must be 
axisymmetric, Commun.Math.Phys. 271 (2007) 

699-722, 'gr-qc/0605106 . 

[34] V. Moncrief and J. Isenberg, Symmetries of Higher 
Dimensional Black Holes, Class. Quant. Grav. 25 (2008) 
195015, arXiv:0805.1451 . 

[35] S. Hollands and A. Ishibashi, Black hole uniqueness 
theorems in higher dimensional spacetimes, 

Class. Quant. Grav. 29 (2012) 163001, 

[arXiv: 1206.1164. 

[36] S. Green, S. Hollands, A. Ishibashi, and R. M. Wald, 

In preparation. 

[37] O. J. Dias, G. T. Horowitz, and J. E. Santos, Black 
holes with only one Killing field, JHEP 1107 (2011) 

115, arXiv:1105.4167 . 

[38] O. J. Dias, G. T. Horowitz, D. Marolf, and J. E. Santos, 
On the Nonlinear Stability of Asymptotically Anti-de 
Sitter Solutions, Class. Quant. Grav. 29 (2012) 235019, 
arXiv:1208.5772 . 

[39] O. J. Dias, J. E. Santos, and K. Skenderis, “Black 
moons.” In preparation. 

[40] also known as the Garter solution. 

[41] The onset of superradiant modes with oppositely signed 
amplitudes yield the same physical black resonators. 

[42] Similar black resonators may also exist inside Dirichlet 
walls with E and J conserving boundary conditions. 




